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We study a PT-symmetric quantum mechanical model with an 0(A^)-symmetric potential of the 
form m?x^ /2 — g(x^Y /N using its equivalent Hermitian form. Although the corresponding classical 
model has finite-energy trajectories that escape to infinity, the spectrum of the quantum theory is 
proven to consist only of bound states for all A'^. We show that the model has two distinct phases in 
the large- A'' limit, with different scaling behaviors as goes to infinity. The two phases are separated 
by a first-order phase transition at a critical value of the dimensionless parameter rr? j g^^^ , given 
by 3-2^/^ 
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I. INTRODUCTION 



Models with PT symmetry have emerged as an interesting extension of conventional quantum mechanics. There is 
a large class of models that are not Hermitian, but nevertheless have real spectra as a consequence of PT-symmetry. 
Bender and Boettcher have shown that single-component quantum mechanical models with PT-symmetric potentials 
of the form —A {—ixf have real spectra [1]. An extensive literature on PT-symmetry and related matters now exists, 
and there are extensive review articles available f2| 1^. Most of the results have been for models without continuous 
internal symmetries, but there have been some results on models with 0{N) symmetry (l) [5] H]. These models are 
particularly interesting in the large- iV limit. Meisinger and Ogilvie have shown that a PT-symmetric version of the 
O (iV)-invariant anharmonic oscillator is isospectral with a Hermitian model with an 0{N — 1) symmetry [5]. We 
study here the properties of this model, using its Hermitian form. 

The Euclidean Lagrangian of the PT-symmetric model with 0{N) symmetry is given by 
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where g is positive. The minus sign in front of g would lead to a Hamiltonian unbounded from below if the model 
were Hermitian. Prom the standpoint of PT symmetry, the interaction term can be considered as a member of a 
family of PT-invariant interactions 
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which are invariant under PT symmetry [l]. This class of models is well-defined for p = 1, and must be defined for 
p > 1 by an appropriate analytic continuation of the Xj as necessary. In [5], it was shown that this PT-symmetric 
model is equivalent to a Hermitian model with Euclidean Lagrangian given by 

7f + -ma + —a -t- —a tt - ^2gNa (3) 
where cr is a single variable, and tt is a vector of — 1 variables. The corresponding Hamiltonian H is 

This Hermitian form has many remarkable features. First of all, it has a manifest 0{N — 1) symmetry associated 
with rotations of tt rather than the 0{N) symmetry of equation (ml. As in the similar case of a PT-symmetric —gx'^ 
theory of a single variable [7] |8], there is a linear anomaly term which breaks the classical symmetry a —a at order 
h. The 7? field has no quadratic mass term, while the sign of the mass term for the a field is opposite the sign of the 
X fields in the original Lagrangian. From a naive field-theoretic point of view, the tt field is massless at tree level. 

There are two important questions we will address. The first question is the presence or absence of scattering 
states in this model for any value of A^. It is clear that the corresponding classical model has a class of trajectories 
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with finite, continuously varying energies that escape to infinity along paths with tr = 0. Nevertheless, we will prove 
that the quantum mechanical model has only bound states with discrete energy levels. This behavior is similar to 
that of the quantum-mechanical x^y^ model. This model can be derived from the quantum-mechanical reduction of 
a two-dimensional gauge theory |9] and was originally thought to be fully ergodic, i.e., to have only chaotic motion. 
This turns out not to be the case [lOj. Simon has shown that the quantum mechanical version of the x^ip model has 
a purely discrete spectrum jTl]. As we will show in section |llj the arguments of Simon can be generalized to show 
that the PT-symmetric 0{N) models have discrete spectrum for all finite values of N. 

The second, more difficult question, concerns the existence and interpretation of the large- iV limit. We explain 
the nature of the problem in section III. In section IV, we show that a simple variational approximation gives us 
the clues we need to prove that this model has two distinct scaling behaviors in the large- limit, controlled by 
the dimensionless parameter rr? jg^^^ . A first-order transition occurs in the large- iV limit when rr? j g^l^ — 3 • 2^1'^ . 
Because this is a quantum mechanical system, the phase transition only appears in the strict limit of iV ^ oo. For 
large but finite A^, there is a rapid crossover between the two different scaling behaviors. Section V rederives the 
results of section IV using the Born-Oppenheimer approximation, which proves to be exact in the large-iV limit. 
In section VI we present the results of a numerical study of the ground state energy for rv? — that confirm our 
analytical results. A final section presents our conclusions. 



II. ABSENCE OF SCATTERING STATES 

We now turn to the issue of the spectrum for finite TV. In [11], Simon gave five arguments for the absence of 
scattering states in an x'^y^ potential model. His first, and simplest, argument is based on a lower bound for the 
Hamiltonian 
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using the operator inequality 



dx'^ 



+ x'^y'^ > \y\ 



(6) 



which obviously holds on a harmonic oscillator basis. Applying this inequality symmetrically to x and y, we see that 
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dx"^ 
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This lower bound Hamiltonian has only bound states. This key step then leads via the Golden-Thompson inequality 
to the conclusion that H^y has purely discrete spectrum. 
Applying the same ideas to our Hamiltonian, we see that 
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The specific bounding Hamiltonian is not so important; the crucial feature is that the bounding Hamiltonian has only 
bound states. This is sufficient to guarantee that H has only bound states for all finite values of N. 



III. SCALING ARGUMENTS AND THE LARGE- iV LIMIT 

We now consider the large- TV limit of the model. Our naive expectation based on the Hermitian 0{N) model is 
that the ground state energy will be proportional to TV as TV 00 [12]. In order to explore this possibility, we rescale 
the Lagrangian by cr ^ y/Na to obtain 

Le = yCT^ + - TVmV^ + 45TVcr'' + Wga^Tr^ - Ny^cr. (9) 

We see that the anomaly term survives in the large- TV limit, unlike the PT-symmetric matrix case [5]. After integrating 
over the TV — 1 tt fields, we have a large- TV effective potential Ve// for a: 

Veff/N = -m^a^ + Aga^ + ^^/32ga^ - ^ga. (10) 
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Figure 1: The effective potential Veff/N versus a for g = 1. Three different values of are shown, corresponding to region 
I, region II, and the critical point. 



The new term ^ ^/^i2go^ comes from the functional determinant for fluctuations of the if fields, and represents the 
zero-point energy of their quantum fluctuations given a constant value for a. The shape of the potential is controlled 
by the dimensionless parameter rv^jg^^^. Figure l] shows the effective potential as a function of a with g set to 
1 for three different values of rr? . We refer to the region where vr? jg^l^ < 3 • 2^/^ as region I, and the region 
where rr? jg^l^ > 3 • 2^/^ as region II. In region I, Veff has a global minimum at cr = 0. For a near 0, Fe// is 
approximately linear, but with different slopes for a > and a < 0. This occurs because the zero-point energy of the 
if fields has virtually the same form as the anomaly term. However, the zero-point energy term respects a discrete, 
classical a — > —a symmetry which the anomaly explicitly breaks. The boundary between regions I and II is given by 
f7^2^^2/3 _ 3 . 2i/3^ where V^ff has two degenerate minima. In region II, V^ff has a global minimum with a ^ 0. This 
change in the behavior of the effective potential as is varied is not seen in the corresponding Hermitian model 
|12) . and naively indicates that the — 1 7? fields are massless modes in region I. 

To understand better the suspect character of the above analysis, it is useful to rederive these results using hy- 
perspherical coordinates |13l [T3] . This formalism will also be used later in the context of the Born-Oppenheimer 

1 /2 

approximation. If we define p — (v?-^) , the reduced, radial Hamiltonian in the sector of angular momentum I is 
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where I is a non-negative integer. After rescaling cr ^ V Na and p v Np, we have 

1 1 92 {N + 21- '3) (N + 21-1) 



rad 



Nm^cr^ + AgNa^ + 16gNa^p^ - 



2N da^ 2N dp^ 8Np^ 
Taking / = 0, we see that the ground state energy to leading order in large- iV is given by minimizing 

1 



TV 



mV^ -I- Aga^ + IGga^p"^ - ^ga 



with respect to p and cr. Minimizing with respect to p, we find 



and must now minimize 



N \-m^a^ + Aga"^ + 2^/2gcj^ - y^a 



(11) 
(12) 

(13) 

(14) 
(15) 



which is identical to our previous expression for V^ff- However, we now notice that when cr = 0, p is infinite. This 
strongly suggests that our treatment of the large- limit is not valid in region I. 



4 



IV. A SIMPLE VARIATIONAL APPROXIMATION 

In order to gain analytical insight about the large- A?^ behavior of this model, we apply a simple variational approx- 
imation using harmonic oscillator ground states. We assume that the ground state has the form 



*[a,7r] = Vo (ct) H ^o{nj) 



(16) 



where cpo is a harmonic oscillator ground state of frequency Q and expected value (tTj) = 0. The wave function ipo (a) 
has frequency w and expected value {a) = v. We then have the variational inequality Eq < Ey^r where 
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V2^v (17) 



provides an upper bound for all v and positive oj and CI. We are free to rescale the variational parameters v, uj, and 
f2 as we like, and we have considered the class of rescalings of the form v — * N^-v, w — > N^ui , and f2 — > N'^Q.. There 
are two different rescalings with non-trivial behavior in the large- TV limit. 

If we rescale the variational parameters as u ^ N^^'^v, to Nuj , leaving f2 unchanged, we find after some algebra 
that the ground state energy in this large- A'' limit is given by minimizing 



E,,„,r = N 
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This in turn reduces to minimizing 



E,, 



N 



(18) 



(19) 



with respect to v. Unsurprisingly, this is equivalent to minimizing our previous expression for the effective potential 
in the conventional large- limit and is valid in region II. 

A different scaling behavior, which we will show is valid in region I in the next section, is obtained if we perform 
rescalings v — > N~^l^v, u) — > N^^^u), and CI — > N~'^/^Cl. This yields a large- A?" limit of the form 



E = N^/^ 
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Minimization of E^ar with respect to the three parameters oj, Ct, and v leads to the solution 
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7Vl/3^1/3 



(20) 
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(22) 
(23) 
(24) 



Numerically, this gives an upper bound on Eq of approximately I.IVJQT N^/^g^/^ . As we show below, this variational 
result is exact in the large- A^ limit in region I. 



THE BORN-OPPENHEIMER APPROXIMATION AND THE LARGE- LIMIT 



We now apply the two different scaling behaviors we have found in the previous section to the Hamiltonian directly. 
If the a field is rescaled by cr ^ N^/'^a, the rescaled Hamiltonian is 



(25) 
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which will turn out to be valid in region II. All of the terms in the Hamiltonian are of ordor A^, except for the 
kinetic energy term for a, which is of order This suggests the use of the Born-Oppenheimer approximation, 

in which heavy degrees of freedom are treated classically. In this approximation, the wave function is written as 
V'// (tt, o) = Uii (tt, a) wii (cr) where un (tt, a) satisfies 
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uii (tt , a) = eii (cr) uu (tt , a) 



(26) 



describing the N — 1 n fields as harmonic oscillators with frequencies determined by a. The total energy to leading 
order in N is again the effective potential Ve//- As we have seen, in region II there is a non-trivial solution for which 
the rescaled ground state energy Eq/N has a finite, negative value as ^ oo . In region I, we find that Eq/N ^ as 
N ^ 00 ; the vacuum expectation value of the original, un-rescaled a also obeys N~^/^ (cr) ^ in this limit as well. 
The other rescaling is cr —> N~^/^(j combined with tt — > N^/^tt. The rescaled Hamiltonian is 



H = ArV3 



1 



1 



(27) 



With this rescaling, it appears energies will increase with N as N^/'^. We also see from the kinetic energy term for 
cr that the cr field retains its operator character as N goes to infinity, unlike the conventional large- limit obtained 
in region II. We can immediately drop the and terms as irrelevant in the large- A' limit. The relevance of the tt 
kinetic field is less clear. If we again use hyperspherical coordinates for tt, we find that the effective Hamiltonian in 
the I = sector can be written for large N as 



H, 



eff 
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(28) 



where p is now the magnitude of the rescaled tt. It is clear that p has a mass of order N"^, and we again use the 
Born-Oppenheimer approximation, this time applied to p. The wave function is written as ipj (tt, cr) = uj (cr, p) wj (p) 
where ui (a, p) satisfies 



'2da 



2 + 16ffa2p2 _ 



ui (ct, p) = ei (p) ui (cr, p) 



(29) 



which describes a particle in a harmonic potential. The ground state energy is 2y/2gp- 1/ (32p2). The energy of the 
combined system is given by minimizing 



1 



32p^ 



+ 2^p 



(30) 



which gives the ground state energy in the large- A?^ limit in region I to be exactly given by 

4/3 



4'^ = (2) A^^V/^ 



(31) 



which is identical to the result of the variational treatment in section IV. 

We now have results for two different large- A?^ limits. In one case, the ground state energy is proportional to A'^ 



E^q"^ = N -min 



J\/32^ 



2ga 



In the other case, the ground state energy 
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7Vl/3ol/3 



(32) 



(33) 



is proportional to N^^^ and is always positive. If Eq^^^ is negative, as it is in region II, it will be favored over Eq'^\ In 

region I, the formula for £^0^^'' appears to predict that the ground state energy is zero. However, this is misleading: it 
is actually predicting that the ground state energy is not growing linearly with A^, but at some less rapid rate. This 
kind of behavior is shown by £^0^'', and the two expressions are in fact consistent. K we examine the behavior of the 
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scaled ground state energy Eq/N in the large- iV limit, it is zero in region I and negative in region II. However, the 
true behavior of Eq in region I is given by Eq^\ Note that similar considerations apply to the expectation value (a). 
In region I, (cr) is decreasing as iV~^/^, while in region II it is increasing as N^^^. 

In the large- limit, these two behaviors are completely incommensurate, and give rise to a first-order transition 
at rr? jg^l^ — 3 • 2^/^, the boundary between regions I and II. However, there can be no phase transition in quantum 
mechanics for any finite number of degrees of freedom. We can understand the behavior of these two solutions for 
finite but large values of N from the effective potential Vejf in Figure [l] In the region of parameter space where 
rr? I g^l'^ is close to 3 • 2^/'^, Vg// has two distinct local minima, and there will be tunneling between these two minima, 
leading to a small splitting of the ground state energy from that of the first excited state. This is very similar to the 
behavior of the double well with a small term linear in the coordinate x added. In this case, however, as TV becomes 
large, the tunneling between the two minima is suppressed, leading to a first-order transition in the large- iV limit. 



VI. THE CASE = 



In order to check our large- iV result in region I, we study numerically the scaling behavior of the ground state 
energy for the case rn} — for large values of N . This value for rr? is in region I, where the ground state energy 
scales as N^l'^ in the large- limit. As we have seen in the previous section, the rr? dependence disappears in this 
limit in region I, so the case rv? = gives the large- behavior throughout this region. After rescaling the fields by 
a g~^^^a and tt g~^^^Tf, the radial Hamiltonian Hrad is given by 



Clearly, g can be set to 1 at — with no loss of generality. Different angular momentum sectors do not mix, and 
we restrict ourselves to the / = sector. 

We have performed a series of variational calculations using a harmonic oscillator basis for both cr and p. In these 
calculations, a finite-dimensional square matrix is formed from the matrix elements of Hmd using a suitable basis. 
The lowest eigenvalue of this matrix is an upper bound on the ground state energy. For p, these basis elements 
were chosen to be solutions of the reduced Hamiltonian for the {N — l)-dimensional harmonic oscillator. With the 

angular frequency scaled to 1, the general form of a basis function for p is p^^'^~^~^'' L^^ri^'^^'^^^p'^) exp [— ^p^] , where 

TO is a non-negative integer, and has energy to -f- {N — 1) /2. Basis sets of size were constructed using n basis 
elements for both a and p . The optimum values for the angular frequencies associated with a and p were determined 
by minimizing the lowest eigenvalue using a nine-element basis. These angular frequencies were then used with a 
36-element basis to estimate the ground state energy. In Figure [2] we plot the 36-element estimate for the ground 
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IGa^p^ 



(34) 
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state energy, scaled by a factor of N , versus N. We also show a fit of the form 

7V-i/^So = a+^. (35) 

Although the fit was obtained using the points from N — 50 to N — 100, the agreement with the calculated energies 
is quite reasonable even at = 5. The constants obtained from this fit are a = 1.71705(4) and b — —2.29(2), where 
the errors are estimated from Richardson extrapolation and varying the fitting conditions. These estimated errors do 
not represent the true systematic errors due, e.g., to truncation to a finite basis. Note the excellent agreement of a 

with the exact result (f)''^^ « 1.71707. 

VII. CONCLUSIONS 

The PT-symmetric 0{N) model has been analyzed in a fairly complete way using its dual Hermitian form. The 
Hermitian form of the Hamiltonian has classical trajectories of finite energy that escape to infinity, and also — 1 
variables which are massless in perturbation theory. Nevertheless, the model can be proven to have only discrete energy 
eigenstates. There are two distinct large- A^ limits, with incommensurate scaling behavior. For rr? jg^l^ > 3 • 2^/^, the 
large- A^ limit is similar in behavior to that of the Hermitian 0{N) model, with a ground state energy depending on 
TO and proportional to A^. For rr? j g^l'^ < 3 • 2^/^, the ground state energy is proportional to A'^/^, and is independent 
of TO throughout this region. The two regions are separated by a first-order phase transition in the limit A^ — > oo. All 
of these results were obtained from the dual Hermitian form of the original PT-symmetric model. Although it was 
shown in [5] that the large- A^ limit of region II could be obtained from the PT-symmetric form, these arguments do 
not appear to give the detailed information for region I that can be obtained from the Hermitian form. It would be 
very desirable to achieve an understanding of this model for both regions using only the PT-symmetric form of the 
model. Such an understanding might be an important step in constructing PT-symmetric scalar field theories. 
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